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We numerically study 2-photon processes using a set of 
harmonics from a TkSapphire laser and in particular interfer- 
ence effects in the Above Threshold Ionization spectra. We 
compare the situation where the harmonic phases are assumed 
locked to the case where they have a random distribution. 
Suggestions for possible experiments, using realistic parame- 
ters are discussed. 

42.50.Dv, 42.50.Ar, 32.80.Rm, 33.80.Rv 



I. INTRODUCTION 

High Order Harmonic generation (HOHG) refers to the 
radiation emitted during the interaction of a short laser 
pulse with a gas jet. The spectrum of such a radiation 
consists of a series of odd harmonics of the fundamen- 
tal frequency having approximately equal conversion ef- 
ficiency (plateau), and followed by an abrupt cut-off. A 
complete description of harmonic generation which al- 
lows for a direct comparison with the experiments should 
include not only the response of the single atom to the 
laser field, but also the propagation of the different har- 
monics generated through the macroscopic medium. Re- 
cently, several studies on harmonic generation have fo- 
cused on the properties of the phase of the generated 
harmonic waves considering both, single atom and prop- 
agation effects . Within the frame of the single atom 
response the strength and the phase of each harmonic 
wave corresponds to the modulus and the phase acquired 
by the dipolc moment d(t), or equivalently, the dipole ac- 
celeration induced by the external laser field. Thus the 
spectrum of the radiated harmonics can be directly cal- 
culated from the Fourier transform of the dipole acceler- 
ation 101: 



D{uj) = J d{t)e lut dt = A(u 



(1) 



where A(oj) is the amplitude, 4>{lu) is the phase of the 
dipole acceleration in the frequency domain, and |Z?(w)| 2 
the strength of the harmonics. The single atom response 
shows that the phase of each harmonic, 4>{uj), is usually 
shifted with respect to the fundamental phase and de- 
pends strongly on the intensity of the laser field 0] . Let 



us recall that the harmonics are said phase-locked if the 
phase of each individual harmonic is a linear function of 
its order q so it can be written as <j) q = (3quj + C where 
(3qoj is the linear q-dependent term and £ a constant term 
identical for all harmonics. In general, harmonics from 
the plateau region present a random phase when com- 
pared to each other, and only the harmonics from the 
cut-off region are locked in phase. This can be intu- 
itively understood using the semi-classical description of 
the atom-field interaction in which harmonic generation 
results from the recombination of those electrons which 
have previously tunneled out through the effective poten- 
tial produced by the Coulomb and the laser field. The 
induced dipole, and hence the process of harmonic gen- 
eration corresponds to trajectories in which the electron 
returns to the core with appropriate kinetic energy and 
recombines emitting a photon. Each contribution con- 
tains then a phase factor equal to the real part of the 
action acquired by the electron following the respective 
trajectory Among all possible trajectories, very few 
of them are relevant. For the harmonics in the cutoff 
there is only one dominant trajectory, corresponding to 
a single recombination time and, therefore, the phase of 
the harmonics in this region remains locked. On the other 
hand, for the harmonics lying in the plateau, there are 
typically two relevant trajectories corresponding to two 
different returning times. The phase of such an harmonic 
contains the contributions of these two trajectories and 
in particular their quantum interferences. This quantum 
intcrcfcncc term leads to an apparently random phase. 
In spite of this fact, when one analyzes the harmonic 
emission in time rather than in frequency domain, one 
observes a clear periodicity of the signal. Antoine et al. 
Q have shown that the time dependent emission consist 
of a train of ultrashort pulses with two dominant pulses 
per half cycle, each corresponding to one of the relevant 
trajectories. Further studies show that under certain ge- 
ometrical conditions it should be possible to phase-match 
the contribution of only one of these trajectories. That 
results in a train of ultrashort pulses equivalent to those 
obtained by combining several plateau harmonics with 
their phases locked. 

From an experimental point of view, determining the 
relative phase of the emitted harmonics is not a straight- 
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forward task. Its study necessarily involves some type 
of interference effects between the different harmonics. 
Although some pioneering work on interference between 
harmonics of the same order have been recently reported 
a conclusive result on the relative phase has not 
yet been achieved. 

In this paper we address this question using a combi- 
nation of successive odd harmonics to study two-photon 
processes. When several harmonics irradiate simultane- 
ously an atom, there exist various paths associated with 
2-photon transitions ending up into the same final state 
thus leading to possible interference effects. Those effects 
should be clearly present in the ionization rates, as well 
as in the photoelectron spectra. Similar effects of sensi- 
tivity to harmonics phases have been previously studied 
in a series of papers on ATI by mixing harmonics and in- 
frared laser fields Furthermore, two-photon ioniza- 
tion of atoms with photons pertaining to the VUV-XUV 
domain have been recently reported for argon using the 
third harmonic of KrF laser |l0[ ], and for helium using 
the 9 th harmonic of a Ti:Sapphire laser giving thus 
evidence of its experimental feasibility. In this paper, we 
focus on two-photon ionization of atoms by Ti:Sapphire 
harmonics whose order varies from q = 11 to q = 33. 

Essentially, the aim of the paper is to show that mea- 
surable quantities such as ionization or the more de- 
tailed photoelectron spectrum have a different behavior 
depending on the relative phases of the harmonics and 
to study these effects quantitatively. The paper is or- 
ganized as follows: in Section |n| we develop a combi- 
natorial model derived from perturbation theory which 
is able to provide rough scaling laws for the total ion- 
ization rates as a function of the number of harmonics 
present in the field. Section III reports on numerical 
Above Threshold Ionization rates obtained by solving the 
Time-Dependent Schrodinger Equation (TDSE) with re- 
alistic polychromatic pulses for atomic H. Section IV is 



harmonics: 



a study of quasi-resonant multiphoton ionization (MPI) 
in He + also obtained by solving TDSE. Finally, in Sec- 
tion IV we give concluding remarks and we comment on 



possible experiments. 



II. A COMBINATORIAL MODEL 

The model presented in this section provides an esti- 
mate for the ionization probability as a function of the 
number of successive odd harmonics included in the laser 
field. It is based on a combinatorial approach and a per- 
turbative treatment of the atom-field interaction. 

Let us first sketch the problem under investigation. 
We consider electrons emitted during ionization via a 2- 
photon process from an atom initially in its ground state 
(the possible case where some excited intermediate state 
come into resonance with the absorption of only one pho- 
ton is discussed in Section [V ) . Since the field considered 
here consists of a linear combination of successive odd 



E(t) = ^ E q cos(qu> L t + <j) q ) 



(2) 



photons of different energies will be involved in the pro- 
duction of electrons with different kinetic energies, as it 
is illustrated on Fig. [j]. Note that the energies at which 
electrons are released appear on the photo-electron spec- 
trum separated exactly by the energy of 2 photons of 
the fundamental field (u>l) which generates the harmon- 
ics. Moreover, electrons at a particular energy may have 
been released by the absorption of photons originated 
from distinct harmonics. The only restriction being that 
whatever the combination of photons is, the sum of their 
energies ought to be the same. For example, the produc- 
tion of electrons in the central peak of the electron spec- 
trum of Fig. |l| results from transitions following three 
different quantum paths. The amplitudes associated to 
each quantum path may interfere depending on the rel- 
ative phase of the photons involved. The interference 
pattern appears clearly if one writes the ionization prob- 
ability associated with these processes (labeled 3 on Fig. 
|l|) within perturbation theory: 

P/b = \ a M i f 2 g ) E 1 E 3 e l ^ + ^ + 
rt 2 )/ 



b Mf g >(E 2 e l,p2 ) 2 



M^EsE^^+^l 2 , 



(3) 



where the subscripts 3,1 and 2 refer to the biggest, the 
smallest and the intermediate photon, E refers to the 
maximum field amplitude, and <f> refers to the phase of 
the photons. The transition amplitude matrix elements 
from the initial state \g) to the final state |/) for the 
three paths are given by: 



a M (2} = lim 
J 9 



(f\z\n)(n\z\g) 



U>1 



M 



(2) = Unl ^ (Mn)(n\z\g) 



(2) 
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(f\z\n) (n\z\g) 
uj g + uj 3 - uj n + ie 



(4) 



where the summation runs over all coupled intermedi- 
ate states including the continuum. The special case of 
ATI can be numerically extrapolated p| or computed 
accounting explicitely for the resonant free-free transi- 
tion matrix elements |13|. Since we assume the interac- 
tion with a linearly polarized field and express it in the 
length gauge, only the matrix elements of z appear in 
expression (0). The ionization probability then reads: 



P fg =\ a Mg\ 2 EfEl 
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ElElEle-^ 4 ' 1 ^-' 1 ^ +c.c). (5) 
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Interferences in the probability given by Eq. (||) are gov- 
erned by the phase difference A = 0i + 03 — 202 that 
is, the relative phase between the harmonics involved in 
this process. Such a partial probability can be derived 
for all accessible final states. Note, that interferences are 
always absent in processes involving less than 3 different 
photons which is always the case for the processes leading 
to the two first and two last final states (like processes 
labeled 1,2,4 and 5 in Fig. |l|) independently of the total 
number of harmonics included in the field. 

When the incident field is a linear combination of suc- 
cessive harmonics, there exist several distinct quantum 
paths associated to transitions to the same final state. 
These paths may interfere depending on the relative 
phase between the photons. Obviously, the number of 
possible interfering paths increases with the number of 
different frequencies included in the field. 

To simplify the problem, we will now assume that all 
harmonics have the same strength, that is Ei = E 2 = 
E% = E. As the harmonics we are considering here all 
come from the plateau part of the spectrum, this assump- 
tion is generally fulfilled. Furthermore, we also assume 
that all 2-photon transition amplitude matrix elements 
(as they are defined in Eq. (g)) are real and have the 
same value, i.e. a M fg = b M fg = c M fg = M . Although 
never fulfilled in general, this approximation is reason- 
able if the first photon absorbed brings the system into a 
region free of intermediate resonances. For the case of our 
illustrative example, the ionization probability associated 
to the central process (labeled 3 in Fig. [l]), corresponds 
to: 



p fg = M 2 E 4 (5 + 4cos(0! + 3 - 2</. 2 )). 



(0) 



Let us now analyze the two configurations of harmonic 
phases we are interested in, that is, the phase locked case 
and the random distribution of phases. As defined in the 
introduction, if harmonics are phase locked to each other, 
their phase has the form: <f> q = flqu + £. Applied to our 
example A = 01+03 — 202 simply reduces to zero. There- 
fore, the cosine function in Eq.(@) takes value one. On 
another hand, we can interpret the random distribution 
of phases as a loss of coherence between different quan- 
tum paths. It is then reasonable to associate the random 
distribution to the average value of Eq.@ over all pos- 
sible phases. Since the average value of the cosine func- 
tion is zero, all terms in the probability vanish except the 
constant (phase-independent) terms. Generalizing these 
results for an arbitrary number TV of incident harmonics, 
we obtain that the ionization probability associated to 
the central peak behaves as: 



P( PL \N) = a 2 N 2 



and 



P, 



(HP) 



(N) 



27V - 1 if TV odd 
2N if TV even 



(7) 



(8) 



where PL refers to phase locked configuration, RP to 
random phase configuration and a = ME 2 . The to- 
tal 2-photon ionization probability is straightforwardly 
obtained by integrating the electron spectrum over the 
electron energy range relative to the 2-photon processes. 
This probability accounts for the side peaks in the photo- 
electron spectrum and is obtained by simply adding the 
probability associated to each side peak (see Fig. [j]). It 
can be shown that the total probabilities read: 



(9) 



and 



, /27V 3 TV 

a 2 1 

3 3 



p( RP HN) = a 2 [4 N{N - 1) +N 



(10) 



a 2 ( 2N 2 - N) 



As mentioned earlier, when ordered by increasing final 
state energy, the two first and two last processes will 
never show path interferences independently of the total 
number of harmonics involved. Also, the process having 
the largest number of interference terms will always be 
the central one (like process number 3 in Fig. [j]). 

To summarize, the total ionization associated with 2- 
photon absorption should vary approximately as 2/3 A 3 if 
the harmonics are locked in phase, and like 27V 2 if the rel- 
ative phases between them is random. Accordingly, the 
ionization rate associated to the central process should 
vary like TV 2 in the phase locked case and like 27V if the 
phase distribution is random. 

We expect that similar scaling properties will charac- 
terize other processes of interest, such as 2-photon Above 
Threshold Ionization (ATI) which we shall discuss in the 
next section. The advantage of studying ATI instead of 
2-photon ionization processes is that such processes are 
free from intermediate resonances that modify the inter- 
ference pattern, as we shall discuss in Section IV. The 



disadvantage on the other hand, is that the ATI signal 
might be low and, consequently the effects due to the rel- 
ative phase of the harmonics might be difficult to observe 
experimentally. 

Finally, it is worth stressing here that the problem con- 
sidered in this section has a lot in common with the prob- 
lems of ionization by stochastic laser fields, discussed in 
the literature in the 70's and 80's [|l4j. In particular, the 
enhancement by a factor Kl of the A-photon ionization 
probability by multimode laser fields has been predicted 
and observed. In our random phase model, the same ef- 
fect is responsible for the fact that p( flp ) ~ 2!7V 2 . On 
the other hand, the problem considered in this paper is 
more general and deals with the differences induced by 
random and phase locked configurations in laser pulses in 
real atoms, i.e. it fully accounts for atomic bound states 
and continuum structures. 
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III. 2-PHOTON ABOVE THRESHOLD 
IONIZATION 

In this section, ATI is investigated using a time- 
dependent non-perturbative approach. It consists 
in solving the Time-Dependent Schrodinger equation 
(TDSE) for a hydrogen atom interacting with the com- 
bination of several harmonic fields. The method used 
to solve the TDSE and to compute the electron spec- 
trum is based on the expansion of the total wave- function 
on spherical harmonics and B-spline functions (see [jig] ). 
The temporal propagation in the length gauge is well 
adapted to the rather low intensities considered here. 
The electron spectra correspond to electrons emitted in 
the direction of the polarization. 

A. Ionization of hydrogen in its ground state 

We have computed the electron spectra obtained by 
shining a linear combination of successive harmonics on 
atomic hydrogen in its ground state. The frequencies 
of the selected harmonics are chosen so that absorption 
of a single photon already brings the system into the 
continuum. Thus, the 2-photon processes we are inter- 
ested in belong to the so-called Above Threshold Ion- 
ization regime. In our approach, the field is described 
semi-classically as: 

N 

E{t) = E f(t) cos{quj L t + <p q ), q = 2(n - 1) + q Q , 

71=1 

(11) 

where f(t) is the normalized envelope of the field, in our 
case a cosine squared with 5 fs full width at half max- 
imum (20 cycles of harmonic 11); E D is the maximum 
field amplitude (identical for all harmonics) correspond- 
ing to an intensity of 10 13 W/cm 2 ; Ulol — 1-5 eV is 
the fundamental photon energy and q indicates the low- 
est harmonic component in the polychromatic field. The 
number of harmonics we consider varies from N = 2 to 
N = 5, so that the odd harmonics included are H13- 
H15 for N = 2, H13-H17 for 2V = 3, H11-H17 for N = 4 
and H11-H19 for N = 5. The field corresponding to the 
case N = 5 is plotted in Fig. ||. The electron spectrum 
due to such a field is shown on Fig. || where we observe 
that electrons produced via different order processes are 
clearly distinguishable. The spectrum has to be read as 
follows: there are 3 clear subsets of ATI peaks, each of 
them corresponding to a different order process. The left 
most subset corresponds to order 1, that is, direct photo- 
ionization via the absorption of a single photon (Hll for 
the very first peak and H19 for the fifth). The second 
subset, which is the one investigated here embodies all 
processes of order 2. For example, the third peak of that 
subset results from the simultaneous absorption of either 



two photons from the harmonic 13, or one from the har- 
monic 11 and another from the harmonic 15. All possi- 
ble allowed polychromatic combinations are included in 
this type of simulation. Higher order processes appear 
with exponentially decreasing probability and, therefore, 
their study have a restricted interest from an experimen- 
tal point of view. 

We have limited ourselves to a maximum of 5 combined 
harmonics in this study because higher values would lead 
to overlapping subsets of ATI peaks as it is almost the 
case for N = 5 where the 3-photon processes subset 
partly overlap the 2-photon subset (see Fig. ||). 

Let us now focus on the 2-photon Above-Threshold 
Ionization, that is, the central subset of ATI peaks in 
Fig. [| The problem of interest here is to determine 
how much does the total ionization associated with 2- 
photon absorption increase as we increase the number of 
harmonics in the field assuming that either all harmonics 
are locked in phase, or their relative phases are random. 

We have first simulated the case of phase-locked har- 
monics. This is achieved by solving the TDSE when all 
harmonic phases <pq are se t 1° zero. Fig. ^ shows the 2- 
photon ionization probability as a function of the number 
of harmonics included in the field. The main result to be 
read from Fig. || is that ionization certainly scales like 
the cubic power of the number of harmonics as predicted 
by the combinatorial approach Eq. ([)]). Similarly, the ion- 
ization probability associated with the central ATI peak 
scales like 2V 2 " 1 , in good agreement with Eq.(Q). 

To simulate the case where the phases are random, 
we have accumulated the results of 500 different random 
phase configurations after solving the Schrodinger equa- 
tion for each of them and we have averaged the results 
at the end. This has been done for each value of N . We 
report in Fig. the ionization probability thus obtained 
together with its statistics fluctuations. The most inter- 
esting result here is the behavior of the average value of 
the ionization (open circles in Fig. ^|). The average value 
scales like N 2A , departing significantly from the combi- 
natorial estimate that predicts a N 2 law (Eq. [l0|) . Let us 
now examine the statistical fluctuations. On one hand, 
we find that the maximum value of the ionization (down 
triangles in Fig. ||) scales as 2V 2 - 9 , very close to that of 
the phase locked case. Obviously, the configurations that 
lead to such maximal value are the special cases where all 
phase differences are close to zero like in our simulation 
of the phase locked configuration. On the other hand, 
the most globally destructive phase configuration lead- 
ing to the minimum value of the ionization probability 
scales like 2V 1 - 9 (up triangles in Fig. ||). We have also 
investigated the ionization probabilities associated to the 
central 2-photon ATI peak, since it is the transition that 
involves the largest number of interfering quantum paths. 
The results are presented in Fig. g. The average value of 
the ionization probability scales as (2V ) in comparison 
with the fit (2V 2 - 1 ) obtained for the phase locked case. 
Finally, its minimum value drops dramatically since in- 
terferences now can be completely destructive, that is, 
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there exists a particular configuration of phases able to 
totally suppress ionization for this particular electron en- 
ergy. This feature is absent in the case of total 2-photon 
ionization since, as explained before, the 4 outermost ATI 
peaks correspond to transitions involving only two quan- 
tum paths, and therefore are not affected by any change 
in phase. 

A way to improve the statistical results is to increase 
the number of harmonics included in the field by using 
higher order harmonics. This will prevent us from hav- 
ing different ATI order processes overlapping but also the 
ionization signal will decrease because the cross-section 
associated to these processes decays with increasing pho- 
ton energy. 

B. Ionization of hydrogen in the 2s state 

A solution to ease 2-photon ionization due to higher 
harmonics is to start from a less bound atomic energy 
level like H(2s) for example. 

By choosing harmonic 21 (Tilol = 1.5 eV) as the central 
component, we can combine up to 7 harmonics. Because 
the lower binding energy, we set now the field intensity of 
all harmonic components to 10 10 W/cm 2 . As expected 
with this intensity, the probability for ionization by a 2- 
photon process is very low. However, one should keep 
in mind that this probability varies like I 2 in agreement 
with the process order, and that the intensity we are 
using here underestimates experimental attainable fluxes. 

We find that in the phase locked case, the scaling laws 
concerning total ionization associated with 2-photon ab- 
sorption shown in Fig. [?] agree very well with the combi- 
natorial estimates expressed in Eqs. (Q) and (j^). 

The random phase simulation for H(2s) has been car- 
ried out by accumulating results from 300 resolutions of 
the TDSE, each time with a new random distribution of 
harmonic phases. In Fig. ||, we show the 2-photon ATI 
probability, obtained by integrating the spectrum over 
the whole range of electron energies corresponding to 2- 
photon ATI processes. In the case N = 7, the integration 
runs from Ek = 40 eV up to Ek = 76 cV. Here again, 
there is a significant departure of the average value of 
the ionization probability, which scales as A 2 - 55 , from 
the N 2 combinatorial law associated with the random 
phase configuration. 

Finally, as in the previous section, we have looked at 
the most phase sensitive ATI peak located at the cen- 
ter of the 2-photon ATI subset in the electron spectrum. 
The expected combinatorial power laws are 2 and 1 for 
harmonic phases, locked and random respectively. The 
quantum mechanical computation gives the values of 2 
and 1.3 for the respective configurations as it is shown in 
Fig. ^. The lower limit of the distribution drops abruptly 
until it reaches the background induced by the numerical 
accuracy (around 10~ 16 on Fig. ||). 



To summarize, we have shown that there is a differ- 
ent behavior in the 2-photon ionization probabilities as 
a function of N if one considers the harmonics locked in 
phase or if the distribution of their phases is random. To 
a rough extend, scaling laws can be explained using the 
idealized combinatorial model of Section II. 



IV. 2-PHOTON IONIZATION 

The previous section has been devoted to the study of 
2-photon ATI processes, where absorption from a single 
photon already brings the system into the continuum. 
Now, we propose to study direct 2-photon harmonic ion- 
ization (MPI). In contrast with the 2-photon ATI cases, 
ionization requires now two harmonic photons, so the 
interference effects governed by the harmonic phases ap- 
pear directly in the first subset of peaks in the electron 
spectra. 

At this point, let us recall that our combinatorial model 
is based on the assumption that the total transition am- 
plitude towards a given final state is only related to the 
number of interfering paths leading to such final state. 
This assumption, in turn, implies that all the matrix el- 
ements involved in the transitions are considered equal. 
The results obtained so far indicate that this is a rea- 
sonable assumption for 2-photon ATI processes where 
bound-bound transitions do not play any essential role. 
However, 2-photon ionization processes (MPI) may in- 
volve quasi-resonant bound-bound transitions. At low 
intensities, i.e. in the frame of perturbation theory, the 
values of the bound-bound transition matrix elements de- 
pend strongly on how close the transition is from a res- 
onance. As we shall see, this fact substantially modifies 
the results we have encountered so far. 

In order to optimize the study of the interferences be- 
tween different harmonics having, at the same time, a 
good control of the role of the resonances involved, we 
choose as a model a hydrogenic ion (He + ). Note, how- 
ever, that the results we present in this section can be 
straightforwardly generalized to any atom or ion and, 
therefore, have a certain generic validity. 

We numerically solve the TDSE for direct 2-photon 
ionization, where the external field is a linear combina- 
tion of different harmonic fields, all of them with the 
same polarization and intensity. As in the previous sec- 
tions, we consider harmonics from a TkSapphire laser 
(JvjJl = 1.5 eV), set the intensity of each harmonic to 
1.2 x 10 12 W/cm 2 and calculate ionization probabilities 
as well as photoelectron spectra for phase locked and ran- 
dom phase configurations. For our choice of parameters, 
the smallest harmonic from Ti:Sapphire required for 2- 
photon detachment of He + is the 19th and the largest 
one the 35th. Larger harmonics will directly detach the 
ion, and smaller ones than 19th cannot achieve 2-photon 
detachment. 

Before proceeding further, it is illustrative to show 
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some of the effects due to atomic resonances in the photo- 
electron spectrum. To this aim we display the spectrum 
for a phase locked configuration containing: (a) N=4 har- 
monics ranging from 19th to 25th (Fig.|lOj) and; (b) N=6 
harmonics ranging from 19th to 29th (Fig.|TT|) . In both 
spectra there are two distinguishable set of peaks: the 
first one corresponds to 2-photon detachment processes, 
the second one corresponds to 3-photon processes (ATI 
peaks). From an experimental point of view, the interest 
in the second subset should be restricted due to the low 
probability of 3-photon processes. 

The first subset contains 2N — 1 different peaks, where 
N is the number of harmonics present in the external 
field. As mentioned in the previous sections, the cen- 
tral peak of this subset corresponds, in both cases, to the 
transition involving the largest number of different quan- 
tum paths. While in Fig. |l0| this peak is the most proba- 
ble one, in Fig. |ll| the central peak is clearly suppressed 
in comparison with lower order peaks. The reason for 
this peak suppression is due to the fact that the 27th 
harmonics (only present in case (b)) crosses the Is — 2p 
resonance quasi-resonantly. 

Let us now focus on the results corresponding to ion- 
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ization probabilities. As in section [II. the 2-photon de- 



tachment probability is calculated by integrating the area 
below all 2-photon processes in the photoelectron spec- 
trum. We start by analyzing the random phase config- 
uration where we have accumulated results from a large 
number of random phase distributions and we have av- 
eraged them at the end. In Fig. [l| we plot the average 
value of the ionization probability versus the number of 
incident harmonics N which varies from TV = 2 to N = 7. 
The lowest harmonic involved is always the 19 th so N = 2 
includes harmonics 19-21, N — 3 includes harmonics 19- 
21-23, and so on. The overall behavior of the ionization 
probability obtained numerically (dots) agrees relatively 
well with the combinatorial approach (full line). 

In Fig. ^ we present the corresponding ionization 
probability as a function of the number of harmonics N, 
for the phase locked configuration. In the simulation, 
all the individual harmonic phases, 4> q , have been set to 
zero. Again, the full line represents the results from the 
combinatorial approach (Eq. |^) , while symbols show our 
numerical results. We observe that for low order inci- 
dent harmonics i.e. from 19th to 25th, the numerical 
results agree relatively well with the combinatorial ap- 
proach. However, when higher harmonics are included 
(from 27 th onwards) the combinatorial approach clearly 
overestimates the ionization probability obtained numer- 
ically. The discrepancy between the estimate and the 
numerical value becomes then dramatic. 

This discrepancy can be understood by analyzing the 
role of resonances within perturbation theory. Under the 
parameters we have used in the simulations, the energy of 
the 27 th harmonic is slightly detuned (from above) from 
the Is — 2p resonance. As a result, the matrix elements 
corresponding to transitions of the type: 



where i refers to the second photon absorbed, flip their 
sign compared to those transitions which are far off res- 
onance. 

Let us analyze in more detail the 2-photon MPI tran- 
sitions. From the 2 possible paths starting from the 
ground state leading to photodetachment; s — > p — > s 
and s — > p — > d, the second one clearly dominates over 
the first one. This is because transitions that increase 
their principal quantum number and their angular mo- 
mentum in the same direction are always favored com- 
pared to the others |l(|. Therefore, we focus on 2-photon 
transitions of the type s — > p — > d and present the values 
of M in table | 

The table has to be interpreted as follows: rows refer 
to the first absorbed photon and columns to the second 
one. For example, the first row displays the values of the 
matrix elements corresponding to transitions in which 
the first absorbed photon is the 19 th harmonic and the 
second one the one indicated by the column label, i.e. 
transition amplitude of the type: Yl ^[j^+J "^[^ EigEi. 

As it should be, we observe that the values of the ma- 
trix elements (i) decrease as the energy of the final state 
increases and (ii) decrease for transitions far from reso- 
nances. Immediately after a resonance has been crossed, 
the matrix element flip its sign. Thus, with the param- 
eters we have use in the simulation, all matrix elements 
corresponding to the absorption from the ground state 
of a photon from the 27 th harmonic (Eq. (|12|)) flip their 
sign compared to all the other matrix elements shown in 
the table. 

This relative sign flip strongly affects the total ion- 
ization probability in the phase locked case, for which 
p(PL) ^ |^Afi| 2 , i.e. proportional to the square of 
the sum of the individual matrix elements. It has also 
consequences for random phase configurations but there 
the effects are less dramatic. It is easier to understand 
the difference between both cases with an example. Con- 
sider an incident field containing harmonics 25 th , 27 th and 
29 th . Such an incident field, leads to five possible final 
electron energy states, which ordered by increasing en- 
ergy correspond to: (i) 25+25 (absorption of 2 photons 
from harmonic 25 th ); (ii) 25 + 27, 27 + 25; (iii) 25 + 29, 
27 + 27, 29 + 25; (iv) 27 + 29, 29 + 27 and finally (v) 
29 + 29. Remember that in order to have interference 
effects we need at least 3 different quantum paths con- 
tributing to the process. Therefore, the partial ioniza- 
tion probabilities associated to states (i), (ii), (iv) and 
(v) are the same regardless the relative phase of the har- 
monics. The partial ionization probability associated to 
state (iii) which contains contributions from 3 different 
quantum paths depends on the relative phases of the har- 
monics. For phase locked harmonics, it is proportional to 
1^(25+29) + ^(27+27) + ^(29+25) ? ■ On the other hand, if 
the harmonics have a random phase, the partial ioniza- 
tion probability associated to the final state (iii) is pro- 
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portionalto |M ( 25+29)+M(29+25)| 2 + |Af 2 7+27| 2 ■ Compar- 
ing both expressions, we note that it is only in the phase 
locked case where the sign of M( 2 7+27) plays a definite 
role since it acts effectively as a destructive interference 
term. Notice that the values of the matrix elements cor- 
responding to transitions 27 + 25 and 27 + 29 do modify 
the partial ionization probabilities associated to the final 
states (ii) and (iv) respectively, but they do it exactly in 
the same way for both cases: locked and random phase 
configurations. 

To corroborate our findings we have extended our nu- 
merical results up to the 33 rd harmonic (N = 8). For 
the parameters we have used, the photon energy of har- 
monic 33 rd is slightly higher than the atomic transition 
Is — > 3p. Consequently, the matrix elements correspond- 
ing to absorption from the ground state of a photon from 
harmonic 33 rd will again flip their sign, and again act 
as an effective destructive interference for locked phase 
configurations. As a result, the disagreement between 
the ionization probability and the combinatorial estimate 
should increase. Our numerical results support this con- 
clusion. 

In spite of the discrepancy between our numerical 
results and the combinatorial estimates, it should be 
stressed here that nevertheless the overall behavior of 
the ionization probability for random and phase locked 
configurations remains distinct, allowing in principle for 
a distinction between both cases. 

We would like to remark here that accordingly to what 
we have seen, it might also occur that for a particular 
final state energy, its associated partial transition prob- 
ability becomes larger for incident harmonics with ran- 
dom phases than for phase locked harmonics. Such a 
partial peak inversion demonstrates that the effects due 
to atomic resonances can be in some cases strong enough 
to compensate the effects due to quantum phase interfer- 
ence. 

To summarize, we have shown that that the overall be- 
havior of the ionization probability clearly depends on the 
relative phase configuration of the incident harmonics. 
Furthermore, we have shown that for phase locked con- 
figurations and incident harmonic fields near-resonance, 
the total ionization probability strongly departs from the 
combinatorial estimate and it does no longer follow the 
iV 3 power law (Eq. ^). If, on the contrary, the incident 
harmonics are all far of atomic resonances, the combina- 
torial estimate predicts accurately the exact results (Eq. 
^ and We stress, however, that most two-photon 

ionization experiments involving a set of harmonics is 
bound to show the effects of resonances. Finally, for in- 
cident harmonics with random phases, the average value 
of the ionization probability agrees relatively well with 
the combinatorial approach i.e.iV 2 . 



V. CONCLUSIONS 

In conclusion, we have shown that in 2-photon ATI 
of hydrogen atoms by multiple orders high harmonics, 
measurable quantities such as the ionization yield or the 
photoelectron spectra significantly depend on the relative 
phase of the harmonics. It should therefore be possible to 
determine experimentally the phase configuration of the 
high harmonics and, in particular if the phases are locked. 
The method suggested by the calculation is a measure- 
ment of these quantities as a function of the number of 
selected harmonics. Clearly this is a difficult task. One 
difficulty lies in the relatively low photon harmonic flux 
currently achievable: 10 9 photons using 100 fs lasers in 
argon, and 10 7 photons using 30 fs in neon Such 
a flux is in general too low to observe two-photon pro- 
cesses. However, several possibilities to increase the har- 
monic fluxes are under investigation, based on shorten- 
ing the duration of the pump pulses fllj,p||ri| or improving 
the phase-matching |p0 |. Furthermore, initial states pre- 
pared as coherent superpositions of two or more bound 
states can potentially increase the harmonic generation 
efficiency in the part of the spectrum we are interested 
in (ljjl^]. Another challenge to be faced in a real ex- 
periment will be the making of multilayer mirrors with 
different bandwidths to vary the number of selected har- 
monics without changing their relative phases. In spite 
of these obstacles, we think that such an experiment is 
not out of reach and could be useful in asserting the exis- 
tence of subfemtosecond harmonic pulses in the cases of 
locked phases. 

Acknowledgments. A.S would like to thank Richard 
Tai'eb for very illuminating ideas concerning section 
IV and acknowledges financial support from DGICYT 
(Spain) contract number PB95-0778-C02-02 and Euro- 
pean TMR program (FMRX-CT96-0080). E.C. thanks 
Pascal Salieres for a careful reading of the manuscript 
and P. Lambropoulos for clarifying discussions. 



[1] For a recent review see P. Salieres, A. L'Huillier, P. An- 
toine, and M. Lewenstein, Adv. At. Mol. Opt. Phys. 1, 1 
(1998), in press. 

[2] K. J. Schafer and K. C. Kulander, Phys. Rev. Lett. 78, 
683 (1997). 

[3] M. Lewenstein, P. Salieres, and A. L.'Huillier, Phys. Rev. 

A 52, 4747 (1995). 
[4] P. Antoine, A. L'Huillier, and M. Lewenstein, Phys. Rev. 

Lett. 77, 1234 (1996). 
[5] R. Zerne, C. Altucci, M. Bellini, M. Gaarde, T. W. 

Hansen, A. L'Huillier, C. Lyng0a, and C. G. Wahlstrom, 

Phys. Rev. Lett. 79, 1006 (1997). 
[6] T. Ditmire and R. Smith, Opt. Lett. 23, 618 (1998). 
[7] M. Bellini, C. Lyng0a, A. Tozzi, M. Gaarde, T. W. 



7 



Hansch, A. L'Huillier, and C. G. Wahlstrom, Phys. Rev. 

Lett. 81, 297 (1998). 
[8] V. Veniard, R. Ta'ieb, and A. Maquet, Phys. Rev. Lett. 

74, 4161 (1995). 
[9] V. Veniard, R. Ta'ieb, and A. Maquet, Phys. Rev. A 54, 

721, (1996). 

[10] D. Xenakis, O. Faucher, D. Charalambidis, and C. Fo- 

takis, J. Phys. B 29, L457 (1996). 
[11] Y. Kobayashi, T. Sekikawa, Y. Nabekawa, and S. Watan- 

abe, Opt. Lett. 23, 64 (1998). 
[12] E. Cormier and P. Lambropoulos, J. Phys. B 28, 5043 

(1995) . 

[13] E. Cormier, H. Bachau and J. Zhang, J. Phys. B 26, 4449 
(1993). 

[14] See for instance: B.W. Shore The Theory of coherent 
atomic excitation, Vol. 2. (J.Wiley & Sons, New York, 
1990). 

[15] E. Cormier and P. Lambropoulos, J. Phys. B 29, 2667 

(1996) . 

[16] H. A. Bethe and E. E. Salpeter, Quantum Mechanics 
of one and two Electron Atoms (Academic, New- York, 
1998). 

[17] Z. H. Chang, A. Rudquist, H. Wang, M. Murnane, and H. 
Kapteyn, Phys. Rev. Lett. 79, 2967 (1997). P. Salieres, 
P. Antoine, and M. Lewenstein, submitted to Phys. Rev. 
Lett. (1998). 

[18] J. Watson, A. Sanpera, X. Chen, and K. Burnett, Phys. 

Rev. A 53, R1962 (1996). 
[19] A. Sanpera, J. Watson, M. Lewenstein, and K. Burnett, 

Phys. Rev. A 54, 4320 (1997). 
[20] A. Rundquist, C. G. Durfee III, Z. Chang, C. Heme, 

M. M. Murnane, and H. C. Kapteyn, Science 280, 1412 

(1998). 



TABLES 







Matrix elements: s — > 


p — > d 








19 


21 


23 


25 


27 


29 


31 


19 


-84.1 


-66.1 


-58.2 


-49.1 


-41.8 


-35.8 


-30.9 


21 


-90.5 


-74.4 


-61.9 


-52.0 


-44.2 


-37.9 


-32.7 


23 


-105.4 


-86.2 


-71.4 


-59.9 


-50.8 


-42.7 


-37.5 


25 


-168.1 


-136.4 


-112.4 


-94.0 


-79.5 


-68.0 


-58.6 


27 


132.1 


105.1 


85.3 


70.5 


59.1 


50.2 


43.1 


29 


-20.6 


-19.9 


-16.9 


-14.5 


-12.5 


-10.9 


-9.5 


31 


-183.0 


-151.7 


-125.3 


-105.2 


-89.5 


-76.8 


-66.6 



TABLE I. 2-photon transition matrix-elements as a func- 
tion of the order of each photon. 
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FIG. 1. Sketch of the processes investigated in the case of 
the mixing of three harmonics. Also shown is the correspond- 
ing electron spectrum. 
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FIG. 2. Electromagnetic field amplitude corresponding to 
a linear combination of 5 odd harmonics ranging from order 
11 to 19 (fiiOL = 1.5 eV). 
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FIG. 3. ATI spectrum of H(ls) ionized by the selection of 
harmonics defined in Fig. ^| The 3 different subsets of ATI 
peaks corresponds to different order of processes. 
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FIG. 4. 2-photon ATI total probability (filled circles) and 
2-photon ATI probability restricted to the central peak (open 
circles) as a function of the number of harmonics N for the 
locked phase case. The line and dashed line are power law 
fitting curves with a respective argument of 3 and 2.1. 
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FIG. 5. 2-photon ATI total probability statistics (500 runs) 
as a function of N. The relative phases are random. Open 
circles refer to the average value of the distribution (power 
law fit with an argument of 2.4. Also shown are the maxi- 
mum, minimum values of the distribution and the standard 
deviation (Triangles down, up and error bars respectively). 
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FIG. 6. Same as Fig. H but restricted to the central peak. 
Open circles: average value. The power law fit gives an argu- 
ment of 0.7. 
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FIG. 7. Ionization from H(2s): 2-photon ATI total proba- 
bility (filled circles) and 2-photon ATI probability restricted 
to the central peak (open circles) as a function of N. The 
relative phases are locked (zero). The line and dashed line 
are power law fitting curves with a respective value of 2.98 
and 2.04. 
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FIG. 9. Same as figure N but restricted to the central peak. 
Open circles: average value. The power law fit gives an argu- 
ment of 1.3. 
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FIG. 8. Ionization from H(2s): 2-photon ATI total proba- 
bility statistics (300 runs) as a function of TV. The relative 
phases are random. Open circles refer to the average value. 
The power law fit gives an argument of 2.5. Also shown are 
the maximum, minimum values and the standard deviation 
(Triangles down, up and error bars). 
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FIG. 10. Photoelectron spectrum of He + for a phase locked 
configuration containing 4 harmonics (from 19th to 25th). 
The first subset corresponds to direct 2-photon ionization 
processes. At higher energies a second subset of peaks corre- 
sponding to 3-photon process appears. 




Photoekdxon energy (a.u.) 

FIG. 11. Same as Fig. hd but for an external field con- 
taining now 6 harmonics components (from 19th to 29th). 
Observe that the peak distribution in the 2-photon ionization 
processes does not follows the same pattern as before. 



FIG. 12. Average value of the ionization probability versus 
the number of harmonic components for random phase config- 
uration. Full line corresponds to the combinatorial estimate 
(Eq.(|lo|)) and dots display the numerical results. 
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FIG. 13. Same as Fig. n2j but for a locked phase config- 
uration. Full line corresponds to the combinatorial estimate 
(Eq.(0)) and dots display the numerical results. 
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